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ABSTRACT 

d We report on the discovery of a fourth component in the HD 221356 star system, previously 

known to be formed by an F8V, slightly metal-poor primary ([Fe/H] — —0.26), and a distant 
I M8V + L3 V pair. In our ongoing common proper motion search based on VISTA Hemisphere 

^ Survey (VHS) and 2MASS catalogues, we have detected a faint / = 13.76 ±0.04 mag) co- 

,__( moving companion of the F8 star located at angular separation of 12. 13 ±0.1 8 arcsec (position 

angle of 221.8±1.7deg), corresponding to a projected distance of ^312au at 26 pc. Near- 
er infrared spectroscopy of the new companion, covering the 1.5-2.4 micron wavelength range 

with a resolving power of R^600, indicates an LI ± 1 spectral type. Using evolutionary models 
the mass of the new companion is estimated at ~0.08 solar masses, which places the object 
close to the stellar-substellar borderline. This multiple system provides an interesting example 
psj of objects with masses slightly above and below the hydrogen burning mass limit. The low 

I mass companions of HD 221356 have slightly bluer colours than field dwarfs with similar 

spectral type, which is likely a consequence of the sub-solar metallicity of the system. 
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1 INTRODUCTION 

Because of progressive cooling with age, brown dwarfs do not obey 
a unique mass-luminosity relation ( [Burrows e t al. 1997| |2001| l. 
Therefore, the determination of a brown dwarf mass requires either 
a good knowledge of its age or a direct dynamical measurement. 
This, in turn, is possible for substellar companions of stars or for 
those that are found in multiple systems. An additional advantage 
is that the metallicity can be inferred from the primary star. For 
solar-type stars the atmospheres are much better understood than 
for very low-mass stars and brown dwarfs, given the poor knowl- 
edge of opacities in cool atmospheres ( [Bean et al.||2006[ |Bonfils| 
[et al. 2005 1. Coeval systems containing low-mass companions also 
provide very useful constraints on evolutionary and atmospheric 
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models (e.g. |Pinfield et al.|2006[|Dupuy et al.|2010| l as well as offer- 
ing a rather unique view on how the process of star formation works 
at the very bottom of the main sequence ( [Burgasser et al.|2007^ . In 
particular, brown dwarf companions with well determined metallic- 
ities are benchmark objects allowing for a better understanding of 
the effects of metallicity on the physical properties and evolution of 
substellar objects ( jPinfield & et al.|2012^ . Unfortunately, substellar 
companions located at wide separations (>50 au) from their parent 
stars are relatively rare, with an estimated frequency of less than 
a few per cent ^ cCarthy & Zu ckerman|2004[ [Lafreniere & et af] 
[20071 [Kraus & Hi llenbrand 2007|^ 

We are conducting a search for very low-mass common proper 
motion companions of nearby (< 25 pc) stars, using the VISTA 
Hemisphere Survey (VHS) (McMahon et al., in preparation) and 
Two Micron All Sky Survey (2MASS) ( [Skrutskie & et al.|2006| l. 
Our sample of objects includes some of the known multiple sys- 
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Table 1. Properties of HD221356A (a.k.a. HR8931, HIP 116106). 



Parameter 


Value 


Refere 


RA. (J2000) 


23''31™31?62 




Dec. (J2000) 


-04°05'16"78 




V (mag) 


6.50 


1 


Spectral type 


F8.0 V 


3 


/;aCos(8) 


178.7 ±0.9 mas/yr 


2 


pis 


- 192.8 ± 0.9 mas/yr 


2 


Parallax 


38.29 ± 0.54 mas 


2 


Distance 


26.12 ±0.37 pc 


2 


Teff 


5976 ± 44 K 


1 


logfe) 


4.31 ± 0.06 cm/s^ 


1 


Mass 


0.94 ± 0.13 M0 


1 


[Fe/H] 


-0.26 ±0.03 


1 


Age 


2.5-7.9Gyr 


1 



l)|Valenti & Fischer ' 2005) 2) van Leeuwen j2007 



3) Caballero(2007l 
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terns ^Faherty et al.|2010| and references therein). One of the targets 
investigated so far was HD 221356, already known to be a triple 
system. The F8.0 V primary is a field star with slightly subso- 
lar metalhcity [Fe/H] = -0.26 ("Valenti & Fischer|2005) , located 
at 26.12±0.37pc (van Lee uwen 2007 j. The main properties of 
this star are given in Table |T] |Gizis et ari ( |2000^ reported that the 
secondary, initially described as a single object, has a photometric 
distance consistent with that of the HD 221356 star determined by 
Hipparcos. The secondary was later resolved by |Close et al.|p002) 
into a binary separated by 0.57arcsec (~14.9 au), using adaptive 
optics on the Gemini North Telescope. Based on their photomet- 
ric colours, they estimated spectral types of M8.0 V and L3.0 V for 
each of the two components, respectively. This binary (hereafter 
referred to as HD221356BC) was also investigated by Caballero' 
( |2007| . Using data at epochs separated by 48.3 years, he confirmed 
a common proper motion between HD221356A and BC. He also 
measured a mean separation of p = 45l"8±0''4 between both 
components, which corresponds to a projected physical separation 
of nearly twelve thousand au, making it one of the widest known 
systems with an L-type component (see also fig 1 1 of |Zhang & et| 
|al.|2010| i. 

In this article we present the identification and characteriza- 
tion of a fourth, very low-mass companion of the HD 221356 sys- 
tem. We outline the procedure and results of our proper motion 
search together with the analysis of / and YJHKs-hanA photometry 
and near-infrared spectroscopic data of the identified companion. 
We derive the physical properties of the new object which turns out 
to be very close to the hydrogen burning mass limit. 



2 IDENTIFICATION AND FOLLOW-UP OBSERVATIONS 
2.1 VISTA Hemisphere Survey data 

The new low-mass companion of HD 221356 A was identified us- 
ing the 2MASS and VHS catalogues. The VHS is a near-infrared 
public survey intended to cover the entire Southern hemisphere 
(--^ 20,000 deg^) in the J Kg broad band filters with a sensitivity 
more than 3 mag deeper than 2MASS. It uses the 4.1-m telescope 
VISTA (Visible and Infrared Survey Telescope for Astronomy) op- 
erating since 2009 at ESO's Cerro Paranal Observatory in Chile 
( [Emerson et al.|2006| l. It is equipped with a wide-field near-infrared 



Figure 1. Proper motion vector-point diagram for the HD 221356 system. 
All correlated objects within 15 arcmin from the primaiy are plotted as 
black dots, with HD 221356 components labelled as A, BC and D. The pri- 
mary is saturated in both surveys, its proper motion value was taken from 
the literature. Time baseline between the 2MASS and VHS epochs is 12.18 
years. 



camera (VIRCAM), comprising sixteen 2kx 2k pixel detectors with 
a mean plate scale of 0.339 arcsec. The HD 221356 system was ob- 
served with VISTA on 2010 November 25 and 26. Average seeing 
conditions were 1.4" and 0.9", respectively. 

The VHS near-infrared images are processed and calibrated 
automatically by a dedicated science pipeline implemented by the 
Cambridge Astronomical Survey Unit (CASU). Standard reduc- 
tion and processing steps include dark and flat-field corrections, 
sky background subtraction, linearity correction, destripe and jitter 
stacking. For more detailed description we refer to the CASU web- 
page http://casu.ast.cam.ac.uk/surveys-projects/vista, as well as to 
[Irwin & et al. ( 2004 ) and Lewis et al. ( 2010 ). 



2.2 Proper motion 

The search for additional common proper motion companions of 
HD 221356 A was done using the astrometry given in the VHS 
and 2MASS catalogues, which provide a 12.18 yr baseline. The 
positions of 2MASS sources have an estimated accuracy of 70-80 
mas over the magnitude range of 9 < A'^ < 16 I Skrutskie & et ali] 
|2b06). The astrometric solution for VHS observations is done au- 
tomatically as part of the CASU pipeline, using the 2MASS point 
source catalogue. The objects on the catalogues extracted from each 
VISTA detector are matched to their counterparts in 2MASS using 
a correlation radius of 1 arcsec. Because 2MASS has a high degree 
of internal consistency it is possible to calibrate the world coordi- 
nate system of VISTA images to better than 0. 1 arcsecond. 

The search was performed using topcatQ a useful tool for 
analysis and manipulation of source catalogues and other data ta- 
bles, developed as part of the Virtual Observatory. We retrieved the 
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Figure 2. False colour VISTA 7-band image of HD 221356AD. Angular 
separation is 12.13 ±0.18 arcsec and the position angle of the identified 
companion is 221 .8 ± 1 .7 degrees. Saturation in the centre of the primary is 
visible. Field-of-view is 1 x 1 arcmin, with North up and East to the left. 

astrometric and photomeric data from both 2MASS and VHS cata- 
logues, for all the objects within a radius of 15 arcmin correspond- 
ing to ~23,000 au around the examined star. To avoid some of the 
spurious detections, we selected sources brighter than 7=17 mag in 
2MASS. We have cross-matched 300 objects from both catalogues 
within 1 arcsec. The sources that remained unmatched were subse- 
quently cross correlated taking into account the proper motion of 
the primary star provided by Hipparcos jvan Leeuwen|2007) . We 
illustrate the resulting proper motion vector-point diagram of HD 
221356 on Figure[T] 

We have found that the proper motion of HD 221356 BC 
(Table |2j is consistent with that of the primary HD 221356 A, 
thereby confirming the result of the Konigstuhl survey of ]Caballero| 
( |2007) . Individual components of the BC pair, separated by only 
0.57 arcsec jClose et al.| 2002 ), were not resolved in the VHS 
images. We also identified a new common proper motion com- 
panion (2MASS J233 13095-0405234, hereafter HD221356D), 
with {^JaC0s8, //g) = (153.48±21,-190.20± 19) mas/yr, located 
12.13±0.18 arcsec southwest from the primary. We adopted a total 
astrometric uncertainty of ~ 28 mas/yr, estimated using the stan- 
dard deviation of proper motions for sources with fj <100 mas/yr. 
The proper motion of the new object is common with that of the HD 
221356 system. Measured separations, position angles and proper 
motions are listed in Table |2] 



2.3 Photometry 

The VHS catalogue provides aperture photometry for HD 221356 
in the YJHKs near-infrared bands. The VISTA photometric sys- 
tem is calibrated using the magnitudes of colour selected 2MASS 
stars converted onto the VISTA system using colour equations, in- 
cluding terms to account for interstellar reddening Photometric 
calibrations are determined to an accuracy of 1-2%. 

The new faint companion HD 221356 D is well resolved in 
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properties 



Table 2. Proper motion, separations and position angles of low mass com- 
ponents of the HD 221356 system. 



Comp. 


fiaCos(5) 
(mas/yr) 


/'S 
(mas/yr) 


p* 
(arcsec) 


8" 
(deg) 


(au) 


BC 
D 


176±21 
154±21 


-167 ±19 
-190 ±19 


451.10±0.18 
12.13±0.18 


261.77±0.04 
221. 8± 1.7 


11900±50 
317±9 


• epoch (MJD) = 55525.12460836; p, £ 


) and r are measured witli respect lo the primary. 



the VHS images (see Figure |2]), but within the glare of the pri- 
mary. In order to minimize the possible light contamination in the 
aperture photometry of the VHS catalogue, we applied a method 
to suppress the PSF of the primary by subtracting the flipped and 
rotated images from the original ones. Standard deviation of the 
background at the same separation as the companion was a factor 
two lower in the PSF subtracted images than in the original ones, 
which slightly improves the detection of the object. We performed 
aperture photometry on the resultant images with an aperture of one 
full-width half-maximum (FWHM). The instrumental magnitudes 
were then calibrated to apparent magnitudes by adding the aper- 
ture corrections determined using the VHS photometry of fourteen 
isolated bright stars located within 10 arcmin from HD 221356A. 
The differences between the new and the catalogue photometry in 
the Y,J,H and K, bands are 0.53 ±0.06, 0.26 ±0.04, 0.15 ±0.03, 
0. 13 ± 0.03 mag, respectively. Photometric values for the BC com- 
panion were taken directly from the VHS source catalogue. 

On November 15th, 2011 we performed follow-up /-band 
imaging of HD 221356 AD. Observations were carried out us- 
ing the IAC80 telescope equipped with a E2V 2048x2048 back 
illuminated CCD detector with a plate scale of 0.304 arcsec/pix, 
which provides a 10.4x 10.4 arcmin field of view. We selected the 
12 images with best seeing (FWHM<1.2 arcsec) to minimize in- 
terference from the primary star. We reduced the images applying 
standard techniques, including bias and flat-field correction, using 
IRAF routines. Individual exposure times were 5s. For each image, 
we performed a similar method for the PSF subtraction of the pri- 
mary as used for VISTA images and subsequently aligned and com- 
bined all of them. We obtained the PSF-fit photometry using the 
DAOPHOT package in IRAF and calibrated the instrumental mag- 
nitude of our object using 1 1 bright stars in the field with DENIS 
l |Epchtein & et al.||l999[ > /-band data available. We note that the 
photometric system used (Cousin) is not the same as the one of DE- 
NIS, and that some differences may appear for very cool objects, 
however in our previous photometric calibrations we found that the 
zeropoint between IAC80 and Denis /-band has small colour de- 
pendence l [Costado et al.|200 5^. 

Additionally, we acquired /-band observations of 
HD 221356 A using FastCam, mounted on the 1.5-m Carlos 
Sanchez Telescope at the Teide Observatory on January 31, 2012. 
FastCam is a lucky imaging instrument, designed to perform high 
spatial and time resolution observations ( |Oscoz & et al.|[2008| l. 
Optics provide a plate scale of 43.5 mas/pix, and a field of view of 
--^ 22 X 22 arcsec^. We obtained 17 blocks of 1000 images of 50 ms 
individual exposure times. Images were bias corrected, aligned 
and stacked into the final image using the software provided by 
the FastCam team. We derived the /-band aperture photometry 
of the primary, since we found that the literature values based on 
photographic plates are not reliable. Instrumental magnitudes were 
calibrated using photometric standard stars from |Landolt| jl992| l 
observed at different airmasses along the night under photometric 
conditions. We also explored the inner region to search for the 
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presence of additional companions, but none were identified. We 
may exclude the presence of an equal mass companion to the 
primary at separations greater than 0.2" (~5 au) and companions 
with A/ < 5 mag at separations greater than 1" (--^26 au). 

The photometric data are listed in Table [3] The /-band mag- 
nitudes of the distant BC pair were taken from jGizis et al.|j2003^ 
who imaged individual components with Hubble Space Telescope 
WFPC2. The integrated JHKs photometry of the BC component 
from VHS were also decomposed into individual magnitudes using 
the flux ratios derived by |Close et al.|p002 . The JHKs band pho- 



tometry of the primary, also given in Table 
catalogue. 



3^ is from the 2MASS 



on a single feature that can be uncertain bearing in mind the poorly 
understood effects of low metallicity. 

We also note that the water vapour absorption band at ~ 1.5 
jAva, in the blue part of the //-band is more intense in this object 
than in early L dwarfs and more resembles mid/late L dwarfs. This 
is not an instrumental effect since it does not appear in the M8.5 
object, but it may be due to a larger contamination by the primary 
at blue wavelengths than at red wavelengths. We therefore cannot, 
with confidence, assign this feature as unusual for this object. In 
summary, we estimate that HD 221356 D is a slightly blue early L 
dwarf with a spectral type of L0-L2. 



2.4 Near-IR spectroscopy 

We also obtained near-infrared spectroscopy of the HD 221356 
AD system and the M8V spectroscopic standard LP 412-31 ( |Kirk-| 
Ipatrick et aT][T995l > using LIRIS, with the HK grism and the 
IKx IK Hawaii detector at the 4.2 m William Herschel Telescope 
(WHT) on December 30 201 1. This instrumental configuration pro- 
vides a nominal dispersion of 9.7 Apix^' and a wavelength cover- 
age of 1.4-2.4 jAm. A slit width of 0.75" was used rotated to the 
direction along the AD system and the final resolution of the spec- 
trum was 26 A(R --^ 600). Total integration time was 2240 s, divided 
into individual exposures of 160s. A nodding pattern of two posi- 
tions (AB) was used to subtract the sky background. Weather con- 
ditions were photometric and the average seeing was 0.9". Data 
were dark corrected, sky subtracted, aligned and combined at each 
nodding position. Flat-field correction using a tungsten lamp was 
not applied due to a spurious feature in the A'-band tungsten spec- 
trum. After subtracting the contribution of the primary wings, spec- 
tra of HD 221356 D were optimally extracted using the APALL rou- 
tine and wavelength calibrated using ArXe arc lines. We finally 
combined the spectra at both AB positions and corrected for tel- 
luric absorption lines by dividing them by the A3V star HR8840, 
observed at a similar airmass, and multiplying by a blackbody 
of the corresponding effective temperature of 8500 K. Spectro- 
scopic data of LP 412-31 were reduced and analized in a similar 
way to HD 221356 D, but telluric correction was done using the 
A3V star HR 1036, also observed at a similar airmass. Spectra of 
HD221356D and LP 412-31 in comparison with other standard 
objects and the main spectroscopic features are shown in Figure[3] 
HD 221356 D displays stronger water vapor absorption bands 
than the standard M8 dwarf observed with the same instrumental 
setup and sky conditions and reduced in the same manner as our 
target. This implies that HD 221356 D has a cooler spectral type, 
quite likely within the L domain. Aimed at deriving the spectral 
type of HD 221356D, we have compared its LIRIS spectrum with 
data extracted from the IRTF libraries ( [Gushing et al.|[2005] l, see 
Fig. 3. We note that these data correspond mostly to dwarfs with so- 
lar composition. The overall HK spectral energy distribution (SED) 
of HD 221356D is better reproduced by a spectral type of LO-Ll. 
However we determine a spectral type of L1-L3 if we consider the 
different water indexes at ~ 1.7 or 2.0 /im ( |Testi & et al.|2001| 
[AUers & et al.|2007||Sle"snick et al.|2004l >. The differences can be 
explained if the object is slightly metal poor (consistent with the 
primary), because the /f-band flux should be reduced by the H2 
coUisional induced absorption. The strength of the sodium feature 
at 2.2 microns in object D compares better with solar metallicity 
M9 dwarf than with the early L-type templates, however we prefer 
to rely on the general SED, which is better fit by LI templates, than 



3 PHYSICAL PROPERTIES OF HD 221356D 

Assuming a distance of 26.12±0.37 pc to the system, we calcu- 
lated the absolute magnitudes of individual components and con- 
structed the {Ml, I — J) and (My, J — K^) colour-magnitude dia- 
grams, shown in Fig. |4] All four objects clearly follow a well de- 
fined photometric sequence, with the new companion located be- 
tween the known B and C components. All the colours and mag- 
nitudes of the new object are in good agreement with its physical 
membership to the HD 221356 system. To better illustrate the po- 
sition of the three low mass companions, we added in both pan- 
els of Fig. |4] the Pleiades low-mass stars and brown dwarfs (at 
d = 120 pc) from Bihain et al. ( 2010 ), M and L dwarfs from Liebert] 
|& Gizis|p006) wit h available parallaxes and the field M dwarfs 
from Leggett et al.| ( |20(j0l l. The Pleiades cluster offers a homoge- 
neous collection of objects with similar age and metallicity. The 
right panel also includes L and T dwarfs from jVrba & et al.] ( |2004[ l, 
[Kirkpatrick et al.| ( |201 1\ , and 12 subdwarfs with spectral types from 
M7 to L7 with measured parallaxes ^Faherty & et al.|20I2| l. In the 
{Mj, I — J) diagram we have also plotted the 5 Gyr isochrones of 
the NextGen models for solar and low metallicity ([M/H] = —0.5) 



Baraffe et al._199 8) and of the DUSTY model for solar metallicity 



Chabrieretal.|2000:l. 



We have depicted a veriety of objects with similar spectral 
types, that span from young ages (--^120 Myr) and solar com- 
position, to old stars with solar metallicities and field subdwarfs 
with metal-poor abundances ([Fe/H]~ —0.5, [Opine et aL]|2007| l. 
The low-mass components of the HD 221356 system have slightly 
bluer colours than typical field stars, placing them on the blue 
edge of the photometric sequence defined by late M and L dwarfs. 
However their colours are redder than known field subdwarfs, im- 
plying an intermediate metallicity, being thus in good agreement 
with the metallicity determination of the primary. The J — and 
I — J colours of the new companion, which are 1.01 ± 0.06 and 
2.94 ±0.14 mag respectively, correspond to the typical values of 
M7-M9 field dwarf s( Leg gett & et a l. 2002, K irkpatrick & Mc-[ 
[Carthy|1994p . These photometric colours suggest an earlier spec- 
tral type than that determined in the spectroscopic analysis. We at- 
tribute this difference to the sub-solar metallicity of the system, 
which affects the spectral energy distribution. In particular, the flux 
suppression in the K^-hzud was already recognized in a group of 
L dwarfs, as a low metallicity feature ( [Kirkpatrick & et al.|2010[ 
[West et al.|2011^ . Because of that, the spectral classification of the 
wide binary components of the HD 221356 system, which is based 
on photometry, may be uncertain. 

The luminosities of the B, C and D components were derived 
from their JHKs-band magnitudes, using the trigonometric distance 
of the primary, bolometric corrections from Golimowski & et al.[ 
( [2004^ and spectral type-colour relations frorri |Vrba & et al.[ ( [2004| l. 
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Table 3. Photometric data for the components of the multiple system HD 221356. 



Component 



/ (mag) 



Y (mag) 



J (mag) 



H (mag) 



Ks (mag) 



A 
BC 
B 
C 
D 



5.953±0.018 
15.536±0.027 
15.568±0.018 
19.363±0.039 
16.70±0.10 



13.695±0.025 



14.933±0.059 



5.488±0.019 
12.852±0.010 
12.933±0.011 
15.713±0.042 
13.763±0.038 



5.264±0.038 
12.261±0.026 
12.353±0.026 
14.993±0.061 
13.209±0.026 



5.150±0.017 
11.946±0.026 
12.055±0.028 
14.495±0.114 
12.755±0.025 



Notes: /-band mag of A is from the FastCam, JHK^ mag are from 2MASS. / mag of the BC i s from|Gizis 
|etal.l(l2003ll, JHK, mags are from VHS and were decomposed using the flux ratios derived by |Close et al, 
I2OO21. / mag of D is from 1AC80 measurement, JHK, mags are from our photometry on VHS images. 
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Figure 3. Near-infrared H- and K,- band spectra of the new companion compared with the M8 standard LP 412-31 jKirkpatrick et al.|1995) , observed with 
the same instrumentation, and the LI (2MASS J14392836-I-1929149) and L3 templates (2MASS J15065441-I-1321060) taken from the IRTF library < Cushing| 
|et al.| 2005 1. Spectra were normalised at 1.7 microns and offsets have been added for clarity. The grey area indicates the region of high telluric absorption. The 
most prominent molecular and atomic features are indicated. 



Table 4. Physical properties of low mass companions in HD 221356 system. 



Comp. 


SpType 


log L/Le 


Teff (K) 


M/Mp 


B 
C 
D 


M8.0±1.5 
L3.0±1.5 
Ll.Oil.O 


-3.28±0.09 
-4.35 ±0.09 
-3.61 ±0.10 


2351-2717 
1756-2141 
2079- 2304 


0.090 ±0.008 

065+° °" 
0.079 ±0.006 



The effective temperature ranges were calculated adopting the tem- 
perature scale for high-gravity field dwarfs given by |Golimowski| 
|& et aL| ( |2004| l, assuming a spectral type range of L0-L2 for the 
new companion and M6.5-M9.5, L1.5-L4.5 for B and C, respec- 
tively. The resulting values are given in Table[4] The masses of the 
new companion and B and C components were estimated using the 
DUSTY model from the Lyon group ( [Chabrier et al.|2000) , which 
is only available for solar metallicity. We adopted a wide range of 
ages, using the 1, 5 and 10 Gyr isochrones. Masses were derived 
from their luminosities by interpolating the mass-luminosity rela- 
tions given in the models. The differences in mass calculated for 
different ages are lower than the errors in mass determination re- 
sulting from the uncertainties of luminosities. For the new com- 
panion we finally adopted a mass of 0.079 ± O.OO6M0, which is 
the average value obtained using JHKs-hand magnitudes, a spec- 
tral type range of LO— L2, and an age range of 1 — 10 Gyr. To take 
account of the differences in models for low metallicity stars, we 
have checked the mass of object D using the 5 Gyr NextGen model 
for [M/H] = -0.5 i jBaraffe et aL]|1998[ l. We obtained a mass of 



0.083 ± 0.002, which is slightly larger, but still within the uncer- 
tainties of that determined for solar metallicity models. 



4 A BENCHMARK SUB-SOLAR METALLICITY 
MULTIPLE SYSTEM 

The age of the primary star in the HD 221356 multiple system was 
estimated by |Valenti & Fischer| l |2005) to be 2.5—7.9 Gyrs based 
on isochrone analysis. The lithium abundance of the primary star 
(log n (Li) = 2.5 in the usual scale of log n{H) = 12) is typical of 
late F-type stars in clusters with ages in the 2—8 Gyr range ( |Sestito| 
|& Randich||20051 ). The chromospheric activity of the star is also 
typical of a moderately old main sequence star ( [Valenti & Fischer] 
[2505| >. 

HD 221356 is a slightly metal-poor stellar system whose com- 
ponents have masses just above and below the hydrogen burning 
limit. For sub-solar metallicities the stellar-brown dwarf border- 
line is expected to be shifted to higher masses, e.g. to --^ 0.079 Mq 
at [M/H] 



-0.5 (0.072 M0 at [M/H] = 0) (Bar affe et al.|1998| l. In 
such a coeval, old system it becomes particularly interesting to in- 
vestigate the lithium abundances of the very low mass components. 
While the M8 star (component B) and the LO— L2 (component D) 
should have fully burnt their original lithium, component C with 
0.065 solar mass may have preserved some amount of the initial 
lithium content. Theoretical models for solar metallicity predict full 
lithium depletion for such a mass, however this may not be the case 
for sub-solar metallicity, since models also predict a less efficient 
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Figure 4. Left panel: M/, 7 — 7 colour-magnitude diagram of the HD 221 356 system. The positions of the four components are marked with points and labelled 
with the corresponding letters. Pleiades low-mass stars and brown dwarfs (crosses) from Biha in et al.H2010) , M and L dwarfs from |Liebeil & GizisH2006) 
with available parallaxes and M dwarfs (dots) from Leggett et al. (20001 are also plotted. The 5 Gyr isochrones of the NextGen models for solar and low 
metallicity stars (Bara ffe et al._1998 l, represented by a dotted and dashed line, respectively, and of the DUSTY model for solar metallicity ^Chabrier et al.| 
|2000^ , shown as a solid line, are also included. The masses in solar mass units from the DUSTY model are also indicated and marked with open circles in the 
corresponding isochrone. Right panel: Mj, J — K of BCD components of the HD 221356 system. We have also added L and T dwarfs from |Vrba & et al.]^2004} 
and |Kirkpatrick et al.H201 \\ and 12 subd warfs with measured parallaxes from |Faherty & et al.H2012) , depicted by squares. The mean L and T near-infrared 
photometric sequence from |Vrba & et al!]j2004^ is represented by a dashed line. 



depletion at low metallicities ( [Chabrier & Baraffe|1997^ . Observa- 
tions of the Li abundance in the three low-mass components of this 
system will constrain both the evolutionary models and the age of 
the system. 

Multi epoch measurements of the system will allow to detect 
the orbital motion of companion D. Although the estimated orbital 
period (assuming circular orbit) is of the order of 5500 years, the 
relative change of the position would be up to '-^ 14 mas/yr, which 
is measurable using modern high spatial resolution imaging (e.g., 
adaptive optics, lucky imaging). The expected semi-amplitude of 
radial velocity variation of the primary induced by the presence of 
companion D will be of the order of 130 m/s (for 90° inclination), 
however the orbital period is too long to allow a full determina- 
tion of the three-dimensional orbit. Maximum annual variations of 
roughly 0.15 m/s are expected, which may be explored with the 
new generation of ultra high precision spectrographs jWilken & et| 
|al.|2012| l. 



5 CONCLUSIONS 

Using the VHS and 2MASS surveys we have identified a new very 
low-mass companion (HD 221356 D) in the slightly metal-poor HD 
221356 system, which thus becomes a quadruple. The new object 
is located at a projected distance of ~ 312 AU from the F8 primary. 
The four components of the system follow a well defined photo- 
metric sequence. From near-infrared spectroscopy we determined 
L0-L2 spectral type for the D companion. Based on evolutionary 
models its mass is estimated at 0.079 ± 0.006 Mq, and its effective 
temperature is in the range from 2100 to 2300 K. The J — and 
I — J colours of the low-mass components are slightly bluer than 



field counterparts of the same spectral type. We interpret this as a 
result of the low metallicity of the system, which may become a 
reference for the spectral classification of metal poor M and L-type 
field objects. Since the distance and metallicity of the HD 221356 
system are well known, the detailed study of its ultracool compan- 
ions, which are located above and below the frontier between stars 
and brown dwarfs, can provide valuable constrains on evolutionary 
models and, in particular, shed light on the properties of objects on 
the transition from stellar to substellar regime. 
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